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We study three subjects on quark confinement in hadrons in SU(3) C lattice QCD. From the 
accurate lattice calculation for more than 300 different patterns of three-quark (3Q) systems, we 
find that the static 3Q potential is well described by Y-Ansatz, i.e., the Coulomb plus Y-type linear 
^Nj ' potential. We also study the multi-quark (4Q, 5Q) potentials in lattice QCD, and find that they are 

» I , well described by the one-gluon-exchange (OGE) Coulomb plus string-theoretical linear potential, 

03 ■ which supports the infrared string picture even for the multi-quarks. The second subject is a lattice- 

QCD determination of the relevant gluonic momentum component for confinement. The string 
tension (confining force) is found to be almost unchanged even after cutting off the high-momentum 
gluon component above 1.5GeV in the Landau gauge. In fact, quark confinement originates from 
£SJ ' the low-momentum gluon below about 1.5GeV. Finally, we consider a possible gauge of QCD for the 

quark potential model, by investigating "instantaneous inter-quark potential" in generalized Landau 
gauge, which describes a continuous change from the Landau gauge to the Coulomb gauge. 

PACS numbers: 12.38.Aw, 12.38.Gc, 12.39.Jh, 12.39.Pn, 14.40. Rt 

43 ; I. INTRODUCTION 

In 1966, Yoichiro Nambu [l[ first proposed the SU(3) C gauge theory, i.e., quantum chromodynamics (QCD), as a 
^ ■ candidate for the fundamental theory of the strong interaction, just after the introduction of "color" Around 
1970, he also proposed the string theory for hadrons 0|, which describes hadron phenomenology. In 1973, the 
asymptotic freedom of QCD was theoretically shown U and QCD was established as the fundamental theory of 
the strong interaction. However, in spite of its simple form, QCD creates thousands of hadrons and leads to various 
• , interesting nonperturbative phenomena such as color confinement and chiral symmetry breaking [j| . Even now, it is 
very difficult to deal with QCD analytically due to its strong-coupling nature at the low-energy. Instead, lattice QCD 
is now a reliable numerical method to analyze nonperturbative QCD. In this paper, we study three subjects on quark 
confinement in hadrons in SU(3) C lattice QCD. 
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II. INFRARED STRING PICTURE FOR B ARYONS /MULTI-QUARKS 

Around 1980, the first application of lattice QCD Monte Carlo simulations Q was done by M. Creutz for the static 
quark-antiquark (QQ) potential. Since then, the study of inter-quark forces has been one of the important issues 
in lattice QCD [7|. Actually, in hadron physics, the inter-quark force can be regarded as an elementary quantity to 
connect the "quark world" to the "hadron world" , and plays an important role to hadron properties. 

Around 2000, we performed the first accurate reliable lattice QCD study for the three-quark (3Q) potential, which 
is responsible to the baryon structure at the quark-gluon level. Furthermore, in 2005, we performed the first lattice 
QCD study for the multi-quark potentials, i.e., 4Q and 5Q potentials, which give essential information for the multi- 
quark hadron physics. Note also that the studies of 3Q and multi-quark potentials are directly related to the quark 
confinement properties in baryons and multi- quark hadrons. 

Here, we summarize our detailed studies of the inter-quark forces in three-quark/ multi-quark systems with SU(3) C 
quenched lattice QCD |Hl0| . 

In lattice QCD, the static QQpotential VQq(r) is calculated with the Wilson loop, and is well described with the 
Coulomb plus linear form as [3, Q 
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This physically means one-dimensional flux-tube formation between quark and antiquark, and this one-dimensional 
squeezing of color flux is shown in lattice QCD 0- 

For the color-singlet baryonic 3Q system, the 3Q potential can be calculated with the 3Q Wilson loop Q. For more 
than 300 different patterns of 3Q systems, we perform accurate calculation of the static 3Q potential V^q in lattice 
QCD on various lattices: (£=5.7, 12 3 x 24), (£=5.8, 16 3 x 32), (£=6.0, 16 3 x 32), (£ = 6.2, 24 4 ). We find that V 3Q 
is well described by the Coulomb plus Y-type linear potential (Y-Ansatz) 



^3Q = -^3Q ^ 1 1 + CT3Q£ m in + C3Q, (2) 

where L m - m is the minimal total length of of the color flux tube linking the quarks, i.e., the Y-shaped flux-tube length in 
most cases. As an example, we show in Fig. 1(a) the 3Q confinement potential V^q nf , i.e., the 3Q potential subtracted 
by the Coulomb part, plotted against the Y-shaped flux-tube length L m i n . At each £, clear linear correspondence 
is found between V^ c q nf and L m ; n , which indicates Y-Ansatz. This physically means Y-shaped flux-tube formation 
among three quarks, and the Y-type flux-tube formation is actually observed in lattice QCD from the measurement 



of the action density in the spatially-fixed 30 system as shown in Fig. 1(b). Y-Ansatz has been also supported 
by many other recent studies in lattice QCD QJfl and theories including AdS / CFT [l3[ . 




FIG. 1: (a) The 3Q confinement potential V 3 C Q nf , i.e., the 3Q potential subtracted by the Coulomb part, plotted against L m i n 
of Y-Ansatz in lattice unit [jj. (b) The lattice QCD result for Y-type flux-tube formation in a spatially- fixed 3Q system [Til ]. 
The distance of each quark from the junction is about 0.5fm. 



After the experimental report of penta-quark candidates, lots of theoretical analyses for the exotic hadrons have 
been done or revisited. Recently, several charmed tetra-quark candidates such as X(3872) and Z + (4430) have been 
experimentally discovered, and tetra-quark has been also investigated as an interesting object in quark-hadron physics. 

Obviously, the quark-model calculation is one of the standard theoretical methods to investigate multi-quark sys- 
tems. In fact, the quark model calculation gives an outline of the properties of multi-quark hadrons and may predict 
new-type exotic hadrons theoretically. However, for such calculations, one needs the quark-model Hamiltonian for 
the multi-quark system. In particular, one has to know the quark confining potential in multi-quarks. Then, we 
performed the first lattice QCD study of multi-quark potential. 

We formulate the multi-quark Wilson loop, and accurately calculate the multi-quark potential in lattice QCD for 
about 200 different multi-quarks [![. We find that the multi-quark potential obeys the OGE Coulomb plus string- 
theoretical linear potential @, 

V nQ = 7^ A nQ j ' 1 + ^nQ^min + C n Q (n = 3, 4, 5, ...), (3) 

z i<j \ Ti v i\ 

where the confinement potential is proportional to the minimal total length L m m of the color flux tube linking the 
quarks. (Figure 2(a) and (b) are examples of the minimal-length flux-tube for the 4Q system.) We find the universality 
of the string tension, o~ n Q ~ 0qq, and the OGE result, A n Q ~ Aqq/2, for n = 3, 4, 5 0, Q. 

For the multi-quark system, the topology of the flux-tube or the string linking quarks can be changed, according to 
the quark location. We also investigate the flux-tube recombination for 4Q systems [|[ . As the example, we show in 
Fig. 2(c) the lattice QCD result of the 4Q potential V4Q for rectangular 4Q configurations as shown in Figs. 2(a) and 
(b), with d = Q1Q2/2 and h = Q1Q3. For large h, the lattice data obey the connected 4Q state energy. For small h, 
the lattice data obey the "two- meson" state energy, i.e., 2Vqq(/i). In fact, the 4Q potential V4Q is found to take the 
smaller energy of the connected 4Q state or the two-meson state In other words, we observe a clear lattice QCD 
evidence of the "flip-flop", i.e., the string recombination between the connected 4Q state and the two-meson state. 
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FIG. 2: (a) A connected 4Q state and (b) a "two-meson" state, (c) An example of the lattice QCD result of the 4Q potential 
V4Q for rectangular 4Q configurations in lattice unit [s|. The symbols denote the lattice QCD data. The solid (dotted-dashed) 
curve denotes connected 4Q (two-meson) state energy. 

In this way, our lattice results physically indicates the validity of infrared string picture, not only for mesons but 
also for baryons and multi-quarks. 

Then, the multi-quark Hamiltonian based on lattice QCD is expressed as 

n n 

^Q = E(P? +M i 2 ) 1/2+ E y OGE(r. i -r, ) + aL min (n = 3, 4, 5, ...), (4) 

i=l i<j 

where Voge is the two-body OGE potential including the spin-dependent part, and Mi the constituent quark mass. 
The confinement part is described by infrared string picture. 

As other recent developments of multi-quark/gluon potentials in lattice QCD, we studied the excited-state 3Q 
potential [l(| and the heavy- heavy-light quark (QQq) potential [lij], and Cardoso-Bicudo studied the three-gluon 
static potential [l5| . 

III. RELEVANT MOMENTUM OF GLUON FOR CONFINEMENT 

Many theoretical physicists consider that confinement is brought by low-energy strong interaction of QCD. However, 
no one has quantitatively showed the relevant energy region for confinement directly from QCD. Here, the key question 
is "What is the important energy scale for confinement in QCD?" Since confinement is mainly brought by gluon 
dynamics, this can be rewritten as "What is the relevant gluon momentum component responsible for confinement?" 

To answer this question, we formulate a new lattice method to extract relevant gluon momentum for each QCD 
phenomenon, based on four-dimensional (4D) discrete Fourier transformation (l6| . For each gauge configuration, we 
select or remove gluon momentum components in lattice QCD by the following five steps. 

1. Generation of link-variable: We generate gauge configurations on a L 4 lattice with the lattice spacing a by the 
lattice-QCD Monte Carlo simulation under space-time periodic boundary conditions, and obtain coordinate- 
space link- variables (x) . We here take the Landau gauge to suppress the artificial gauge fluctuation [l6|, H3] ■ 

2. Fourier transformation: By the 4D discrete Fourier transformation, we define the "momentum-space link- 
variable" , Ufiip) = jj J2 X Uij,(x) exp(z^ i/ p t/ x, y ). The "momentum-space lattice spacing" is a p = 2ir/ (La). 

3. Cut in the momentum space: We introduce a cut (UV cut or IR cut) in the momentum space. Outside the cut, 
Ufi(p) is replaced by the free-field link- variable. We then obtain the momentum-space link- variable U^(jp) with 
a cut. 

4. Inverse Fourier transformation: To return to coordinate space, we carry out the inverse Fourier transforma- 
tion. Since this U'(x) is not an SU(3) matrix, we project it onto an SU(3) element U^(x) by maximizing 
RcTr{J7 A (x) t J7^(x)}. Then, we get coordinate-space link-variable Ufi(x) with the cut. 

5. Calculation: Using the cut link- variable U^(x) instead of U^(x), we calculate physical quantities as the expec- 
tation values in the same way as ordinary lattice calculations. This procedure is applicable to all the quantities 
in lattice QCD. 
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With this method, wc quantitatively determine the relevant gluon momentum component for confinement directly 
from QCD, through the analyses of the QQ potential [l^]. We perform quenched SU(3) C lattice QCD calculations on 
16 4 lattice at 5=5.7, 5.8 and 6.0. 
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FIG. 3: (a) The QQ potential V(R) with the IR cut Air plotted against the inter-quark distance R. (b) The QQ potential with 
the UV cut Auv- The lattice QCD calculation is performed on 16 4 lattice with f3 = 6.0, i.e., a ~ O.lOfm and a p = 2ir/(La) ~ 
0.77GeV. The broken line is the original QQ potential. 

Figure 3 (a) and (b) show the QQ potential V(R) with the IR cutoff Air and the UV cutoff Ayv: respectively fl6j ). 
Then, wc obtain the following latticc-QCD results about the role of gluon momentum components on the inter-quark 
potential. 

• By the IR cutoff Air, as shown in Fig. 3(a), the Coulomb potential seems to be unchanged, but the confinement 
potential is largely reduced. 

• By the UV cutoff Auv, as shown in Fig. 3(b), the Coulomb potential is largely reduced, but the confinement 
potential is almost unchanged. 



p=6.0 1— e— 1 i 
p=5.8 1 — » — 1 

* J 

\%-\ 







0.5 1 1.5 2 2.5 3 3.5 4 



Auv [GeV] 



FIG. 4: The Auv-dependence of the string tension cr(Auv) obtained from the asymptotic slope of the QQ potential V(R) with 
the UV cutoff Auv- The lattice QCD calculations are done on 16 4 lattice with /3 =5.7, 5.8 and 6.0. The vertical error-bar is 
the statistical error, and the horizontal one is the range from the discrete momentum. The broken line denotes the original 
value of the string tension a ~ 0.89GeV/fm. 

Figure 4 shows the Auv-dependence of the string tension ct(Atjv) obtained from the asymptotic slope of the QQ 
potential V(R) with the UV cutoff Auv in the Landau gauge [lg. (Similar results are obtained also in the Coulomb 
gauge.) As a remarkable fact, the string tension (confining force) is almost unchanged even after cutting off the 
high- momentum gluon component above 1.5GeV, while the string tension is significantly reduced when the UV cutoff 
is smaller than about 1.5GcV. We thus conclude that quark confinement originates from the low-momentum gluon 
component below about 1.5 GeV [l6j . 



IV. ATTEMPT OF LINKAGE FROM QCD TO QUARK MODEL 



While QCD is the fundamental gauge theory of strong interaction, the quark model is also a successful model in 
describing hadrons. However, their relation is still unclear. We here consider two paths from QCD to the quark 
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model, based on latticc-QCD results. 



A. Large Gluonic-Excitation Energy and Success of the Quark Model 



In the quark model, most low-lying hadrons can be well described only with quark degrees of freedom. Why 
dynamical gluons do not appear in low-lying hadrons? 

We think that the absence of dynamical gluons is due to the large gluonic-excitation energy of about lGcV for 
mesons [3 andbaryons (See Fig.5.) This is much larger than the quark-origin excitation such as spin-dependent 
interaction. Then, for low-lying hadrons, the gluonic excitation is absent, and the system can be expressed only with 
quark degrees of freedom, which gives a background of the success of the quark model. [In terms of the gluonic 
excitation, the hybrid hadron appears as a higher-excited state of about lGeV, which explains the mass of Y(3940) 
as 4GeV ~ 1.5GeV 2 + lGeV.l 
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FIG. 5: The 1st excited-state 3Q potential Vg 
in lattice QCD at /3 = 5.8 O- A£ 3Q = V 3 e Q S - 



3 e Q ' (filled circle) and the ground-state 3Q potential V 3 g Q S ' (open circle) v.s. L n 
— V^q' is the gluonic-excitation energy. 



B. What is the Gauge of QCD for the Quark Potential Model? 



The quark potential model is a nonrelativistic model with a potential instantaneously acting among quarks. In this 
model, there are no dynamical gluons, and gluonic effects indirectly appear as the instantaneous inter-quark potential. 
From the viewpoint of "gauge" in QCD, the quark model without dynamical gluons can be regarded as an effective 
theory of gauge-fixed QCD. Since the quark model has rotational and global color symmetries, such a gauge should 
keep them like Landau and Coulomb gauges. 




FIG. 6: (a) "Instantaneous potential" V\(R) in generalized Landau gauge for typical values of A. The symbols are the lattice 
QCD results, and each curve denotes the fit with the Coulomb plus linear form, (b) Comparison of the instantaneous potential 
V\(R) at A = 0.02(~ Ac) (black dots) with the physical static inter-quark potential (line). We add the lattice data in Coulomb 
(A = 0) and Landau (A = 1) gauges. 

With a real parameter A > 0, we define generalized Landau gauge (A-gauge) [l9| of diAi + A94A4 = 0. This gauge 
can connect continuously between the Landau (A=l) and the Coulomb (A=0) gauges. In lattice QCD, we investigate 
"instantaneous potential" Vx£R) = -± ln(Tr[Z7|(R, t)U 4 (0, t)]) in A- gauge [l9[ . Figure 6 shows the lattice result of 
the instantaneous potential [19( : no linear part appears in the Landau gauge (A=l); the Coulomb gauge (A=0) shows 
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overconfining with 2^3 times larger confining force [20( , so that the ground state is conjectured to be the gluon-chain 
state [2lj . In the A^-gauge with Xc — 0.02, the physical static inter-quark potential is approximately reproduced by 
the instantaneous potential [19 ]. which seems to match the quark potential model. 
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